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Abstract. We investigate the relation between thermalization following a quantum 
quench and many-body localization in quasiparticle space in terms of the long- 
time full distribution function of physical observables. In particular, expanding on 
our recent work [E. Canovi et al, Phys. Rev. B 83, 094431 (2011)], we focus 
on the long-time behavior of an integrable XXZ chain subject to an integrability- 
breaking perturbation. After a characterization of the breaking of integrability and 
the associated localization/delocalization transition using the level spacing statistics 
and the properties of the eigenstates, we study the effect of integrability-breaking on 
the asymptotic state after a quantum quench of the anisotropy parameter, looking 
at the behavior of the full probability distribution of the transverse and longitudinal 
magnetization of a subsystem. We compare the resulting distributions with those 
obtained in equilibrium at an effective temperature set by the initial energy. We find 
that, while the long time distribution functions appear to always agree qualitatively 
with the equilibrium ones, quantitative agreement is obtained only when integrability 
is fully broken and the relevant eigenstates are diffusive in quasi-particle space. 



PACS numbers: 75.10.Jm, 72.15.Rn, 05.45.Mt 



1. Introduction 



The physics of thermalization in isolated quantum systems has a long and debated 
history. Recent groundbreaking experiments on the non-equilibrium dynamics of low- 
dimensional condensates [TJ [2j [3] triggered a great deal of attention on this topic, which, 
up to them, was mostly addressed as an academic question, in connection with the notion 
of quantum chaos |H El O [TJ El Ej . The demonstration of the lack of thermalization in two 
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colliding bosonic clouds confined in a cigar-shaped potential [3J, and the attribution of 
this observation to quantum integrability, generated a lot of theoretical activity devoted 
to the study of the connections between integrability, ergodicity and thermalization in 
strongly correlated quantum systems [TUl [TT] . The main focus efforts has been on the 
characterization of thermalization resulting from the simplest possible non-equilibrium 
protocol: an abrupt change in time of one Hamiltonian control parameter, that is a 
quantum quench. 

At long times after the quench, the lack of thermalization in an integrable system 
can be seen as a consequence of the sensitivity to the specifics of the initial state, that 
are encoded in the values of the constants of motion along the whole time evolution. 
Following the prescriptions of Jaynes [12], this qualitative understanding was made 
rigorous by the proposal of describing the steady state after a quench by means of a 
generalized Gibbs ensemble, which keeps track of the initial value of all the non trivial 
constants of motion [13]. The conditions of applicability and drawbacks of this approach 
have been extensively tested (see Ref. [10] and references therein). If the system is 
in turn far enough from the integrable limit, thermalization is expected in general to 
occur. This expectation is based on the eigenstate thermalization hypothesis, stating 
that expectation values of few-body observables in a given eigenstate equal thermal 
averages with the corresponding mean energy [U [5] , and it has been tested by means 
of several numerical techniques (see Refs. [TU [TO] and references therein). The issue is 
still under debate [T5| [T6| ITT] . 

A natural scenario to describe the effects of integrability and its breaking on 
thermalization is that of many-body localization. Building on a series of seminal 
papers in disordered electron systems [HI HU [20] , the interplay between integrability 
breaking and many-body localization has been recently studied in the context of 
thermalization [2Tj [22] |23| [21], In analogy to a construction originally conceived 
for disordered electron systems, the quasi-particle space can be thought of as a 
multidimensional lattice where each point is identified by the occupations of the various 
quasi-particle modes. As long as states are localized in quasi-particle space, the 
system behaves as integrable: any initial condition spreads into few sites maintaining 
strong memory of the initial state. Thermalization will not occur. At the same 
time, qualitative behavior of local and non-local operators in the quasi-particle is 
naturally going to be different: locality in quasi-particle space implies sensitivity to 
the localization/delocalization of states, while non-local operators display always an 
effective asymptotic thermal behavior. Once a strong enough integrability-breaking 
perturbation hybridizing the various states is applied, the consequent derealization in 
quasi-particle space will lead to thermalization. An initial state is allowed to diffuse 
into all states in a micro-canonical energy shell generating a cascade of all possible 
lower energy excitations. 

While this scenario appears physically sound, it does not give information on what 
is the degree of sensitivity of the various physical quantities of a many-body quantum 
system to integrability breaking and thermalization. This question is particularly 
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important in view of the fact that recent studies on the dynamics of quantum field 
theories lead to the proposal of a dynamics of thermalization comprising two stages [25] : 
first the system decays to a so-called pre-thermalized state, where the expectation value 
of certain macroscopic observables is to a good approximation "thermal", while the 
distribution function of the elementary degrees of freedom is not [251 ES [27]. At 
a second stage, when the energy is efficiently redistributed by scattering processes, 
real thermalization eventually occurs. Pre-thermalization has been shown to occur 
theoretically for quantum quenches in a variety of systems (26J, ETJ EH EHJ [301 ED 132]). 
Moreover, the study of pre-thermalization in weakly perturbed integrable systems 
shed light on the nature of the pre-thermalized state, which is nothing but a close 
relative of the non-thermal steady state attained asymptotically by its integrable 
counterpart [29]. Signatures of pre-thermalization have been observed in split one 
dimensional condensates [33], which have been shown to be characterized by an 
intermediate, pre-thermalized stationary state. The latter has been investigated by 
studying the full probability distribution of the interference contrast [331 EI] which turns 
out to have a closely thermal behavior, even though the distribution of quasi-particles 
is non-thermal. 

The purpose of this paper is to provide a detailed characterization of the effect 
of integrability breaking on the asymptotic state after a quantum quench, not only by 
studying the average expectation values of certain selected observables, but also focusing 
on their full probability distribution function (FPDF) , following the suggestions that were 
recently put forward in Ref. [M]. We stress that the latter quantity is experimentally 
accessible by studying shot-to-shot variations of a physically measurable observable, 
as it has been performed for the FPDF of matter-wave interference in a coherently 
split Bose gas [33]. In the first part we introduce our model and discuss the long time 
limit of two kind of observables, thus summarizing the results contained in Ref. [22] ; 
in the second part we take a considerable step forward, by extending the discussion to 
quantities which are closer to actual experiments [33], and clarifying more precisely to 
what extent thermalization takes place. We will show that, in the presence of many- 
body localization, the entire probability distribution function describes a canonical 
distribution of the degrees of freedom. 

In the specific, we consider a quantum XXZ spin- 1/2 chain undergoing a sudden 
quench of the anisotropy, in the presence of an integrability breaking term in the form of 
a random transverse field. As the strength of the integrability breaking term is cranked 
up, the many-body level statistics has a well defined transition from Poisson (Integrable) 
to Wigner-Dyson (non-Integrable), closely associated to the localized/diffusive character 
of eigenstates in quasi-particle space [22]. Focusing on the asymptotic state attained 
after a quench from the antiferromagnetic to the critical phase, we compute the FPDF 
of both transverse and longitudinal magnetization densities in a given spatial interval 
by means of exact diagonalization techniques. We compare the results with the FPDF 
of the same observables obtained in a canonical ensemble, with an effective temperature 
fixed by the initial energy. We show that, while for both weak and strong integrability 
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breaking the FPDFs attained after a quench agree qualitatively with those obtained 
in the corresponding canonical ensemble, a full quantitative agreement for the whole 
distribution is only obtained whenever the level statistics in the bulk of the spectrum is 
of Wigner-Dyson type, corresponding to diffusive eigenstates in quasi-particle space. 

This paper is organized as follows. In Sec. |2]we define the models and the details 
of the quantum quench protocol. In the following two sections we first discuss the 
spectrum of the Hamiltonian, which provides an insight in the localization/delocalization 
transition (Sec. [3]), and then we briefly discuss how the long-time asymptotics of different 
correlation functions is affected by the localization/delocalization in many-body space 
(Sec. Hj), thus summarizing the results of Ref. [22]. In Sec. |5jwe take a step further: we 
define the FPDF of the transverse and longitudinal spin and discuss their behavior in 
the asymptotic long-time state attained after a quench (diagonal ensemble). Finally, in 
Sec. [6] we draw our conclusions. 



2. Model 



Throughout this paper we will consider a quantum quench described by a time- 
dependent Hamiltonian: 

U(t)=Ho[g(t)]+'H ib , (1) 

where: 

M J 9o for t < , 
9 ® = [g for t>0. (2) 

The Hamiltonian %{t) is composed of an integrable part %o[g{t)], and an integrability- 
breaking term given by Concerning the integrable part, we will consider the 

anisotropic spin- 1/2 Heisenberg chain of length L (also called the XXZ model) with 
open boundary conditions: 

W z ) =£[ J {afaf ¥1 + afaf +1 ) + J z a^t +1 ] , (3) 

8=1 

where af (a = x,y,z) are the spin-1/2 Pauli matrices on site i, J is the planar xy- 
coupling, while J z is the nearest-neighbor anisotropy parameter in the z direction, 
which coincides with the parameter g that will be quenched at time t = 0. In what 
follows we take H = ks = 1, we adopt J = 1 as the energy scale and work in the 
zero total magnetization sector along the z axis. This model is integrable by Bethe 
Ansatz [35] . The zero-temperature phase diagram is characterized by three regions: a 
gapped ferromagnetic phase ( J z < —1), a gapped antiferromagnetic phase (J z > 1), and 
a gapless critical phase for — 1 < J z < 1. In this gapless region the critical exponents 
depend on J z and the system is characterized by a quasi-long-range order in the xy 
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plane [SB]- I n t ne following we break the integrability of the model by applying a 
random magnetic field: 

L 



where the quantities hi are randomly chosen in the interval [—1, 1]. 

We point out that the model described in Eqs. [3)141 is equivalent to a system of hard- 
core bosons, as one can show by applying the Jordan- Wigner transformation. Therefore 
it is also interesting from an experimental perspective, since cold-atom gases in one 
dimension and at low densities behave like impenetrable bosons [37] . 

After investigating the transition from integrability to non-integrability, we will 
address the long-time behavior of the system following a sudden quench of J z . In this 
context, the knowledge of a substantial part of the spectrum is required, making it 
necessary to resort to a standard exact diagonalization technique. We will diagonalize 
Hamiltonian systems with up to 14 sites, and only consider the zero magnetization 
sector, thus working in Hilbert spaces with up to 3432 states. 

The zero-temperature phase-diagram of the XXZ model in presence of disorder is 
well established [38]. Much less is known at infinite temperature: the phase-diagram 
has been conjectured to be composed of two phases, a non-ergodic many-body localized 
phase (in real space) at A > A cnt , and an ergodic one at A < A cnt , where A cnt ~ 6 4- 8 
at J z — 1 ESj • O ur results indicate the presence of a second non-ergodic localized 
phase (in quasi-particle space) for A = A* close to zero that crosses over to the ergodic 
phase upon increasing A. The fate of this crossover in the thermodynamic limit and 
the eventual value of the critical A* are yet to be determined 0. 

As the strength of A is varied, the system deviates from integrability, as discussed 
more in detail in Sec. [3], where the spectrum and the properties of the eigenstates 
are quantitatively investigated. Indeed, as the level statistics changes from Poissonian 
(typical of integrable systems) to Wigner-Dyson (characterized by the level repulsion 
of non-integrable systems), the eigenstates are also modified. More specifically, when 
the system is integrable they are localized in quasi-particle space, while they become 
diffusive as the system becomes non-integrable. This transition affects the dynamics 
of the system. As summarized in Sec. HI there is a connection between the onset 
of thermalization and the many-body localization transition of the eigenstates (for 
further details, we refer the reader to Ref. [22] ). Later we will go one step forward 
and investigate how the localization/delocalization transition emerges in the FPDF of 
operators [251 ESI E^. 




(4) 



i=l 



| While for the parameters used in this paper the low-lying eigenstates are localized in the 
thermodynamic limit, in the following we consider systems sizes smaller than the localization length. 
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3. Spectral characterization of the integrability-breaking crossover 

In this section we follow the approach of Ref. [22] and consider the properties of the 
eigenvalues and eigenvectors of the Hamiltonian in Eq. [1] for a fixed value of the 
anisotropy J z = 0.5, and show that the integrability breaking is associated with a 
localization/delocalization transition in quasi-particle space. 



3.1. Statistics of the energy level spacings 

In finite-size systems, it is commonly believed that the statistical distribution of the 
energy spacings of the quantum energy levels directly reflects the integrability properties 
of the model [40]. In particular an integrable quantum system is typically signaled 
by the presence of a Poissonian statistics in the distribution of its level spacings {s n }, 
s n = E n+ \ — E n being the spacing between two adjacent levels normalized to the average 
level spacing, 

Pp{s) = e~ s . (5) 

Physically this means that the eigenvalues of the Hamiltonian within a given symmetry 
sector are allowed to cluster. On the contrary, for non-integrable systems level crossing 
is inhibited. The level statistics has a Wigner-Dyson (WD) distribution, which embeds 
the level repulsion in a power-law: lim^o -Pwd(s) ~ s 7 - More precisely, when anti- 
unitary symmetry is preserved, as in the present case, the statistics is described by a 
Gaussian Orthogonal Ensemble [4T] : 



„ . . TVS 7TS^ „ . 

Pwd(s) = —e~- . (6) 

By tuning the parameter A, the XXZ chain undergoes a transition from Poissonian 
to WD statistics. In a finite-size system, this transition takes the form of a smooth 
crossover which can be studied within a specific energy shell, by means of the following 
level spacing indicator (LSI) [8]: 

/ rp\ _ Iq°[P[e,e+w](s) - P P (s)]ds 
Jo [-FWd(s) - Pp{s)\ds 

where P[e,e+w]( s ) is the level statistics computed in the window [E, E + W], while s 
is the first intersection point of Pp(s) and Pwd(s). 

In Fig. [T]we show the level spacing distribution function P(s) for the levels with 
an excitation energy less than a given cutoff , i.e. E < E c , (see caption of Fig. [Q. We 
see that for small A the distribution is closer to an exponential, while for A = 1 it 
almost coincides with a WD distribution, and it is shifted towards larger values of the 
spacings. In the inset of Fig. [1] we show r] w for different intensities of the integrability- 
breaking perturbation. As the strength of the integrability-breaking term increases, the 
LSI approaches values close to unity for A ~ 1. At large values of A ^> J z the system 
tends to another integrable limit [121 H3], indeed we can see that r\ w decreases again for 
A > 1. We note that only the states in the middle of the band display level repulsion 
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(i.e. values of r\ w closer to unity), while states in extreme regions of the spectrum do 
not. This comes as a consequence of the two-body interaction of the system, as opposed 
to the behavior that is typically observed in full random matrices [H]. 

3.2. Inverse participation ratio vs. number of available states 

Along with the Poisson-to-WD transition of the level spacing statistics, the eigenvectors 
also undergo a transition in their statistical properties. Indeed, when the level statistics 
is Poissonian and A is small, the eigenstates of the Hamiltonian are close to those of the 
integrable XXZ chain. In other words, they are localized in quasi-particle space. On 
the contrary, when the statistics is WD, the eigenstates are delocalized in quasi-particle 
space. This picture can be made quantitative by using the inverse participation ratio 
(IPR) [13 SHI S3 SHI HH [50] . Given a pure state and an arbitrary basis {|n)} with 
N elements, the IPR is defined by [41] : [§| 

§ With this definition we fix the typo in Ref. [22] where there was a wrong factor 1 /N in the definition. 
We remark that the results shown in Ref. [22] are not affected by that mistake. 




Figure 1. Main panel: level spacing statistics for three values of A, computed for 
levels with excitation energy E with respect to the ground state up to the cutoff 
energy E c = 20 in our units. The black dotted and dashed-dotted lines show the 
Poisson (Eq. [5]) and the Wigner-Dyson (Eq. [5]) distributions respectively. Inset: LSI 
as defined in Eq. [71 evaluated in a microcanonical shell of width W = 2. Following 
standard techniques adopted in quantum chaos |41j , we performed an unfolding of the 
energy spectrum for each instance. For both plots data are for L = 14 and averages 
are performed over 5000 disorder instances. 
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Figure 2. IPR in the intcgrable basis £j at A = 0.1 (left panel), and at A = 1 (right 
panel), as compared to the number of available eigenstates N in an energy window of 
width W = 2A. Data are for a chain of L = 14 sites. Average over ~ 10 2 disorder 
realizations. 

If a state is a superposition of n st basis states, the corresponding contribution to £ is of 
order n st . Interesting results emerge if one considers two different bases: (i) the "site 
basis" or "computational basis" of the eigenvectors of of with the constraint of zero total 
magnetization, and (ii) the "integrable basis" of the eigenstates of the integrable model 
under investigation (the XXZ chain with J z = 0.5 and A = 0). Given an eigenstate of 
H(J Z ), we found that the IPR in the integrable basis is of the order of unity if A <C J z , 
and it grows with increasing A. Eventually, for large values of A, the eigenstates of 
'H(Jz) become diffusive superpositions, with random phases and amplitudes, of the 
eigenstates of the integrable model. On the contrary, in the site basis, we found an 
opposite behavior, since the states of this basis approach the eigenstates of the system 
at A > J z [22]. 

We can draw an intuitive picture of the localization/delocalization transition 
in quasi-particle space by comparing the IPR in the integrable basis, in a given 
microcanonical shell, with the number of available eigenstates N[e,e+w] m the same 
shell. This is shown in Fig. |2j The microcanonical shell has a width W = 2A ~ V, 
where V is the typical matrix element of the integrablity-breaking perturbation. In 
quasi-integrable situations (A <^ 1, left panel) the IPR is much lower than the available 
microcanonical states, thus meaning that the degree of delocalzation of the system is 
very low. On the contrary, in a chaotic situation (A ~ 1, right panel), the perturbation 
is able to hybridize nearly all the states in the microcanonical energy shell. We point 
out that our approach has been also recently adopted in order to identify the emergence 
of chaos in a very similar spin chain model [50j EI] • 
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Figure 3. Number of states in each microcanonical shell as compared to the initial 
energy of the system. The energies in the x axis has been rescaled on the number 
L of sites, while the number of states N(E) is rescaled with the total dimension of 
the Hilbert space in the zero magnetization sector d(L). Vertical lines denote the 
energy Eq(L) after the quench from an antiferromagnetic ground state, obtained with 
J z q = 100, to J z = 0.5 (the values of L in these vertical lines increase from left to 
right). Different colors refer to different system sizes. We used A = 1 and W = 2A. 
Averages have been performed over ~ 10 3 random instances of the magnetic field. 

4. Long-time dynamics and observables 

4-1. Energy scales involved in the quench 

Let us now consider a sudden quench of the parameter g = J z , from J z q at time 
t < to J z at time t > 0, as in Eq. El We assume that the initial state \ipo) of 
the system is the ground state of %{J z q). Since after the quench the Hamiltonian is 
time-independent, the energy is conserved and is given by Eq = (i/)o\'H(Jz)\i)o) ■ For 
large values of J z q the ground state of the Hamiltonian is the classical antiferromagnetic 
Neel state and the energy density Eq/L converges to a constant value, below the middle 
of the spectral band of the final Hamiltonian. This means that the eigenstates of "H(J 2 ) 
involved in the time evolution of the system are only those lying in the first half of the 
band. This is shown in Fig. [3j where the energy Eq after the quench from an initial 
antiferromagnetic state (vertical lines) are compared with the number of available states 
in a given microcanonical energy shell. 

Following Rossini et al. [521 [53], we then proceed to define an effective temperature 
associated to the quench as the solution of the equality: 

Eq = (H(J z )) TcS = TT[p(T cS )n(J z )] , (9) 
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where p(T cff ) is the equilibrium density matrix at temperature T eff : 

e -H(J z )/T cH 

Equation [9] can be solved numerically for each realization of disorder and then averaged. 
As shown in Ref. [22], the effective temperature increases with the quench strength 
\J Z — J z q\i and tends to saturate for large values of J z , since the ground state approaches 
the antiferromagnetic Neel state. As an example, for a system of L = 12 sites and 
a quench from J z0 = 10 to J z = 0.5 (that are the typical parameters we considered 
hereafter) the effective temperature at A = 1.0 is T e s = 3.4 ± 0.4 (units of J/ks), after 
averaging over 200 disorder instances. 



4-2. Asymptotics of observables 

The connection between the localization/delocalization transition and the onset of 
thermalization deeply affects the dynamics following the quantum quench. A possible 
way to test this interplay is to compare the long-time behavior of the system with the 
thermal behavior at the effective temperature T e g. This can be done at different levels. 
The first possibility is to take a traditional point of view and study the correlation 
functions of a given observable, as we illustrate now. A second more general (to be 
discussed in the following) approach consists in investigating directly the full probability 
distribution function of selected observables. 

Let us start by considering the two-spin correlation functions constructed as 
expectation values of the following operators: 

1 L 

e^V-WLofaf, (a = x, z) . (11) 

j,i=i 

The average value predicted by the canonical ensemble at temperature T e g is given by: 
n^Jk) = (nt) Tee = Tt[p{T eS )nl\. (12) 

The asymptotic value after the quench is found from the diagonal ensemble [HU El |55] : 
n a Q (k) = \\m(mK\m) =X}M 2 M n *l&> , (13) 

^ t— i-OO Z ' 

i 

where \ip(t)) = e~ M ^ J ^ t \ifj ) is the state of the system at time i, while q = (ipol&i) 
is the 2-th component of the initial state |^o) in the basis of the eigenstates 
of the final Hamiltonian %{J Z ). In Fig. |4j we compare the expectation values in the 
diagonal ensemble (black data) of these operators with those predicted by the canonical 
ensemble at the temperature T e g (red data). We have chosen A = 0.4, so that the system 
is still close to integrability and one can appreciate the discrepancies with the thermal 
behavior only for the correlator of cr| (right panel). On the other hand for cr| (left panel) 
a quantitative agreement with the thermal ensemble is observed. Note that the different 
behavior of the two observables is most visible at k = it, where the system is less sensitive 
to boundary effects. As discussed in Ref. [22], the different behaviors of these correlators 
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Figure 4. Comparison between the diagonal (black dots) and canonical (red squares) 
expectation value of the of the two-spin correlation function n x (k) (left panel) and 
n z (k) (right panel) as a function of the momentum k. Data are for L = 12 and 
disorder intensity A = 0.4. Here and in the remaining figures we will always perform 
a quench from J z q — 10 to J z = 0.5. 



at long times could be related to the different nature of the operators involved in the 
low energy limit [52j |53]. Indeed, in the low energy limit [35], where the XXZ chain 
critical phase maps onto a Luttinger liquid and quasi-particles are approximately free 
bosons, cr| turns out to be a local operator, coupling a finite number of quasi-particle 
states, while a% couples all the states and is thus non-local. 

To be more quantitative, we define the absolute "residual" of the operator n£ 
between the diagonal and the canonical ensemble: 

K = I#)-^W|. (14) 

In Fig. El we show such quantity at k = 7r, as a function of A. The residual for the 
non-local operator, 5n*, does not depend substantially on the strength of A nor on the 
size of the system. On the contrary, the residual of the local operator, 5n^, decreases 
significantly as the system departs from integrability and shows a minimum for A = 1. 
For larger values of A the residual 8n* starts to grow. This is consistent with the fact 
that for large values of A the system approaches another integrable limit. Moreover, 
the quantity 8n*, and in particular the position of the minimum A are size-dependent. 
Given the numerical limitations of exact diagonalization, we cannot guess what is the 
limit of A(L) when L — > oo. The value A of the perturbation strength at which the 
system is closer to a thermal behavior, at a given size, corresponds to the situation in 
which derealization in Fock space is most pronounced. 

5. Full probability distribution functions 

The study of FPDFs of observables, instead of only their average value, has attracted a 
lot of interest, both theoretically [3U [56j EH EE] and, very recently, experimentally [33J. 




Figure 5. Residuals 5n% (black curves) and 8n* (red curves) between the diagonal 
and the canonical ensemble predictions. Different symbols refer to different system 
sizes, as depicted in the caption. Averages over 200 instances are performed. 



Here we want to study the probability distribution function of the transverse and of the 
longitudinal spin of a subsystem S of size R < L. 

The xy transverse-spin magnitude of a subsystem S made by R spins is given by: 

2 



(S 



QX 

\°R 



iSl 



° 3 °t) 



(15) 



where Sf 



jes J x j<ies 
^crf is a spin operator on site i and direction a, while spins inside S are 
chosen contiguously in the central portion of the total system and are identified by the 
indexes j = ^ — f+l, . . . , f +f ■ This is exactly the same quantity considered in Ref. [53] . 
where the full distribution of quantum noise in Ramsey interference experiments is 
studied in detail. 

The z component of the spin for a subsystem of size R is given by: 

E^i ( 16 ) 



en 



jes 



Since the eigenvectors of a z R coincide with the states of the computational basis, the 
computational cost of evaluating the probability distribution of this quantity is less 
than that for (S^) 2 . 



5.1. Asymptotics of the FPDFs 



Similarly to what we have done in Sec. I4.2l for the correlation functions, we now consider 
the behavior of the FPDFs in the diagonal [HU [HJ and in the canonical ensemble at the 



Many-body localization and thermalization in the FPDF of observables 



13 



effective temperature T eff , and compare the results. 

The FPDF can be evaluated in the canonical ensemble at the effective temperature 
T e g according to 

P T eS M = E ( Tr K Tcff ) V ln] ) 8{8R - 
n 

= E ^T cff ( S i?,n) ^(Si? - Sfl,n) , (17) 
n 

where PT s { s R,n) = Tr [p(T e s) V^ n ] denotes the thermal expectation value of the 
projector V% n on a given eigenstate of cr£> (a =_L, 2) corresponding to the eigenvalue 
SR, n (see the correspondence with Eq. [T2J). On the other side, analogously to what we 
did in Eq. [131 we define the FPDF in the diagonal ensemble as 

P Q( S R) = E ( E (^1 l&> ) - 

n i 

= E pq( s «,«) - • ( l8 ) 

n 

The comparison between the two ensembles can be made quantitative by computing 
the absolute difference (see EqJT4l): 

5p a (sR,n) = \p a Q (sR,n)-p%JsR,n)\ (19) 

The latter quantity depends on R and on the values of sr specific to each R. We may be 
interested in comparing the FPDF's at different distances: in this case a useful quantity 
is the integrated difference, which is given by: 

Ap a (R) = Y,Sp a (sR,n)/v(R) (20) 

s B.,n 

where Yl' SR denotes the sum over the eigenvalues of compatible with the zero total 
magnetization condition, and v(R) is their number. In the following, we will average all 
these quantities over disorder, and identify their uncertainty with the standard deviation. 
We are going to show our results for a quench from J z0 = 10 to J z = 0.5 in a system 
with L = 12 spins. 

In Fig.|6l[7]and[8]we show the probability distribution of the transverse spin (S^) 2 at 
A = 0.1, A = 0.5 and A = 1.0, respectively. In all the cases, the distribution functions 
have qualitatively the same pattern in both the diagonal and the canonical ensemble. 
But when disorder is small, e.g. for A = 0.1 in Fig. HI there is a visible discrepancy. 
Increasing the intensity of the disorder a quantitative agreement is established, until for 
A = 1, that is, where the system is furthest from integrability, the two distributions 
perfectly overlap (Fig.EJ). The discrepancy between the two distributions is quantified by 
means of the integrated difference Ap ± defined in Eq. [201 which also allows to compare 
the distributions for different values of R. The results are shown in Fig. [9j The common 
feature for all values of R is that, when the disorder intensity is small, i.e. A J z , 
the discrepancy Ap 1 - is larger and more sensitive to the size of R than for bigger A. 
However it is not clear from our results what the dependence on R should look like in 
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Figure 6. Probability distribution of the transverse spin p ± (sR yn ) for a system of 
L = 12 sites and A = 0.1. Different panels refer to different sizes of the subsystem. 
Here and in the remaining figures averages are performed over 500 disorder realizations. 
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Figure 7. Same as Fig. [HI but for A = 0.5. 
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Figure 8. Same as Fig. El but for A = 1.0. 



the thermodynamic limit. Unfortunately with our simulations we can only work with 
very small sizes L, such that the smallest R is only one order of magnitude less that L. 

In Fig. [10] we show the FPDF of the longitudinal spin o~ z R in the diagonal and the 
canonical ensemble for a system of L = 12 sites and R = 6. The constraint Y2j a j = ® 
allows only some eigenvalues of o~ R , so when the size of the subsystem increases from 
R = L/2 to R = L progressively less values of sr are permitted. For this reason we 
choose R = 6, which is the case with the maximum number of eigenvalues allowed by 
the symmetries. We observe in all the cases that there is a maximum at sr = 0, which 
is the most sensitive point to the variation of A. For odd R (data not shown) there are 
instead two symmetrical eigenvalues, sr = ±1, where the function has two maxima and 
which are most suitable to check the dependence against A. From Fig. [10] we see that 
at small values of A there is a difference at sr = between the two ensembles, while 
the FPDFs almost coincide for A ~ 1. The scenario is similar to that of {S R ) 2 . 

In Fig. [TJJwe show the integrated difference Ap z . We see a slight dependence on 
whether R is even or odd. Namely, for even values of R the discrepancy between the 
ensembles at small values of A is more pronounced than for odd R. This effect is more 
visible for small values of A. For even values of R the difference Ap z has its smallest 
value for A = 1, while for odd values of R the minimum is at A ~ 0.5. Looking at the 
error bars, we see that both p z and Ap z are generally characterized by larger fluctuations 
with respect to the transverse spin, perhaps because the spectrum of a R has a higher 
degeneracy with respect to {S R ) 2 . However the error bars, both for the transverse and 




longitudinal spin, are proportional to the disorder intensity. As we already pointed out 
for the transverse spin, we cannot fully characterize the dependence of Ap z on R, because 
we have too small systems. Intuitively, what we can expect is that the constraint on 
the total magnetization £V a j = ® should relate the situations in which the subsystem 
has R and L — R sites respectively. On a qualitative level this is also suggested by the 
panels of Fig. CCD 

6. Conclusions 

In this paper we analyzed the connection between thermalization and many-body 
localization. We studied the long time behavior of a XXZ spin-1/2 chain following a 
quantum quench, when integrability is broken by a random magnetic field. In particular 
we addressed the issue of pre-thermalization, not only looking at the behavior of average 
values of observables, but also at their FPDF. We found a qualitative agreement between 
the FPDF in the asymptotic state and that predicted by a thermal distribution in the 
whole range of disorder intensity we considered. Nevertheless, the situation in which the 
two FPDFs are quantitatively indistinguishable occurs only when the system is furthest 
from integrability, that is when the eigenstates are diffusive superpositions in quasi- 
particle space. 
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Figure 10. Probability distribution function of a z R in the diagonal (black circles) and 
canonical (red squares) ensemble. Data for L = 12 sites, R — 6 and several values of 
the disorder intensity. 




Figure 11. Integrated difference Ap z (R) for a system of L = 12 sites. 
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